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ABSTRACT 

We measure stellar masses and structural parameters for 5,500 quiescent and 20,000 star-forming 
galaxies at 0.3 < z < 1.5 in the Newfirm Medium Band Survey COSMOS and UKIDSS UDS fields. 
We combine these measurements to infer velocity dispersions and determine how the number density 
of galaxies at fixed inferred dispersion, or the Velocity Dispersion Function (VDF), evolves with time 
for each population. We show that the number of galaxies with high velocity dispersions appears to be 
surprisingly stable with time, regardless of their star formation history. Furthermore, the overall VDF 
for star-forming galaxies is constant with redshift, extending down to the lowest velocity dispersions 
probed by this study. The only galaxy population showing strong evolution are quiescent galaxies 
with low inferred dispersions, whose number density increases by a factor of 4 since z — 1.5. This 
build-up leads to an evolution in the quiescent fraction of galaxies such that the threshold dispersion 
above which quiescent galaxies dominate the counts moves to lower velocity dispersion with time. We 
show that our results are qualitatively consistent with a simple model in which star-forming galaxies 
quench and are added to the quiescent population. In order to compensate for the migration into the 
quiescent population, the velocity dispersions of star-forming galaxies must increase, with a rate that 
increases with dispersion. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
fundamental parameters 



1. INTRODUCTION 

The study of galaxy evolution is often focused on 
identifying the primary properties which drive growth 
and eventually determine the properties of any given 
galaxy. Initial studies emphasized the im portance of 
galaxy luminosity (e.g. IBower et al.l 119921 ). Many re- 
cent studies highlight the importance of stellar mass in 
determining both the s tar formation history and the 
black hole growth (e.g. iBundy et al.l 120061 : iPeng et all 
l2010f ). However, hints that stellar mass cannot be 
the sole driver of galaxy evolution remain - for exam- 
ple, at intermediate masses both red and blue galax- 
ies exist (e.g. [S trateva et al.l 120011 : iBlanton et aD 120031 : 
iKauffmann et a l. 2003). Increasing evidence points to 
stellar velocity dispersion, or some combination of stel- 
lar mass and size, as a more fundamental property of 
a galaxy; it is connected to the galaxy's stellar popula- 
tion, the dark matter halo in which it resides and the 
super-massive black hole (SMBH) at its center. Velocity 
dispersion is most notably included as a key parame- 
ter in the Fundamental P l ane for elliptica, ! galax ies (e.g. 
iDiorgovski fc Davis! Il987t iBernardi et al.l 120031) and in 
the iMagorrian et al.l ()1998l ) relation between abuige and 
M,. Recently, striking correlations between velocity dis- 
persion and other galaxy properties such as color and 
star formati on rate s (SFRs) have been examined (e.g. 
Franx et al.1 [20081 : iBell et al.i [20121: Ivan Dokkum et atl 



201 It). FurthermoreT TWake et al.l (|2012aD demonstrated 



using clustering analysis of the SDSS that the stellar ve- 
locity dispersion of a massive galaxy is more strongly 
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related to the halo properties than even stellar or dy- 
namical mass. 

In this context, it is important to create a census 
of galaxies as a function of velocity dispersion and ex- 
plore how it has evolved over time, as it is directly 
linked to the evolution of many facets of galaxy prop- 
erties. The number density of galaxies at a given ve- 
locity dispersion, or the Velocity Dispersion Function 
(VDF) has been studied extensively in the SDSS (e.g . 
Sheth et al."2003: 'Mitchell et al."2005': 'C hoi et al.l [20071: 
Shankar et al. 2009; Bcrnardi et al. 2010) both for ellipti- 
cal galaxies, for which velocity dispersion measurements 
are most commonly used, as w ell as other subgroups 
and the entire galaxy population. iBezanson et al.l (|2011[ ) 
used inferred velocity dispersions to demonstrate that the 
evolution of the VDF out to z ^ 1.5 is differential - with a 
build-up of low dispersion galaxies and little evolution in 
the number of galaxies with high velocity dispersions in 
time, which is qualitatively consistent wit h the expected 
evolu tion of the VDF due to galaxy ages (jShankar et al.l 
120091 ) but so mewhat diff erent from results based on lens- 
ing surveys (|Chad l20lO^. 

Here we extend this work by examining the number 
densities of star-forming and quiescent galaxies sepa- 
rately, with two motivations. First, for the overall pop- 
ulation of galaxies, velocity dispersion (or compactness) 
correlates with quiescence, such that velocit y dispersion 
is anti-correlated with star-formation rat es (jFranx et al.l 
12008*: 'Bell et al.' 2012t IWake et al.ll2012b[ ). Therefore we 
expect the distribution as a function of velocity disper- 
sion to be quite different for each type of galaxy. And sec- 
ond, for individual galaxies, velocity dispersion might be 
a relatively stable property, as opposed to other proper- 
ties such as luminosity, mass and even sizes. However ve- 
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locity dispersion must be able to grow with time, at least 
initially. Based on simple Virial arguments, increasing 
velocity dispersion requires increasing the mass within a 
fixed radiu s and therefore som e sort of dissipative pro- 
cesses (e.g. lHopkins et ani2009[ ). which would mainly im- 
pact star-forming galaxies. On the other hand, it has 
been suggested that velocity dispersions could decrease 
via a process such as minor merging which could puff 
up galaxies without adding a significan t amount of mass 
(e.g. iBezanson et al.ll200alNaab etHI 120091 : lOser et al.l 
mm . Although minor merging could be influential for 
all types of galaxies, this would be more important in 
denser regions or galaxies embedded in more massive 
dark matter halos. 

We present the data used in this analysis in §2 and in- 
vestigate the evolution in the number density of galaxies 
as a function of velocity dispersion and surface density 
in §3. §4 describes a simple model which reproduces 
the qualitative evolution of the VDF. Finally we dis- 
cuss these results and some of the implications for galaxy 
formation in the context of dark matter haloes and the 
growth of SMBHs in §5. Throughout this work, we as- 
sume Ho = 70km/sMpc-\f7M = 0.3 kn^ = 0.7. 



2. DATA 



2.1. 



Redshifts, Stellar Masses, Sizes and Inferred 
Velocity Dispersions 



For this work we use photometric measurements of 
galaxy stellar mass, effective radii and Sersic index to 
estimate inferred v elocity dispersions for galaxies in the 
NMBS CO SMOS (IWhitaker et all [20ll and UKIDSS 
UDS fields (jWilliams et al.ll2009D and study the redshift 
evolution of the VDF. Inferred velocity dispersion has 
been shown to correlate with measured velocity disper- 
sions in the SDSS, especially w hen in corporating the 
gala xy structure in the form of thelSersid ([l 968) n param- 
eter (|Tavlor et al.ll2010l:lBeianson et al.ll2011.) . Details of 
the technique, including the cali bration of this r elatio n 
using the SDSS, can be found in IBezanson et ahl ()2011[ ). 
For each galaxy, we include a variety of optical and 
near-IR broadband and medium b and photometry to fit 
photometric redshifts using EAzYdl rammer "eraLl[2nn8h 
and calculate ste llar masses using FAST dKriek et al.l 
[2001 . assuming IBruzual fc Chariot (2003) stellar pop- 
ulation synthesis models and a Chabricr (2003) Initial 
Mass Function (IMF) and scaling the implied M/L ra- 
tio to the total luminosity of the best-fit Sersic profile. 
By insuring that the luminosity of the galaxy refiects 
the best-fit profile, and not the aperture photometry, 
we can more accurately correct for non-homology such 
that stellar mass better predicts dynamical mass (and 
thereby inferred velocity di spersion will bette r predict 
central ve locity dispersions ) (jTavlor et al.ll20ldl ). We use 
GALFIT (|Peng et al.ll2"nnl to measure sizes using imag- 
ing in two waveband s (UKIDDS J and K for the UDS 
(jWilhams et al.| [200^ and ACS I and WIRDs K for the 
NMBS Cosmos fie ld (Scov ille et aH [2007t IBielbv et al.l 
1201 It IBezanson et al. 2011) ). Utilizing a Sersic depen- 
dent Virial constant (Bertin et al.ll200T 
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and a scale factor fit to SDSS data, we calculate the in- 
ferred velocity dispersion within r^/S in both wavebands 
for all galaxies in the survey. 
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For this analysis, we include galaxies with inferred ve- 
locity dispersions above lOOkm/s, above which the rela- 
tion is calibrated locally and our samples are complete. 
At all redshifts, we compare the inferred velocity disper- 
sions calculated in each bandpass, which agree with a 
scatter of < 0.1 dex at all redshifts. In order to evaluate 
the importance of including the Sersic dependent con- 
stant and correction to total luminosity we also compare 
inferred velocity dispersioi is calculated in each ban dpass 
with k = 5.0 as used in iCappellari et al.l ()2009[ ). We 
find that in addition to increasing the scatter between 
the inferred dispersions in the different bandpasses by 
up to 0.05 dex, eliminating corrections for non-homology 
introduces a net offset such that the K band inferred dis- 
persions are larger by ~ 0.05 dex than those measured 
in the I band ACS imaging. This offset is also apparent 
in the UDS, though the offset is smaller at ^ 0.01 dex 
between the J and K bands. We note, however, that the 
Sersic indices derived from ground-based data (Cosmos 
K band sizes and all sizes in the UDS) may have signifi- 
cant uncertainties (see §5). 

Finally, to limit the effects of color gradients on our re- 
sults, we interpolate the inferred velocity dispersion and 
sizes to a rest-frame wavelength of 6200 A, the central 
wavelength of the r filter where the SDSS inferred dis- 
persions are calibrated. 

2.2. The Effects of Gas 

At the core of this analysis is the fact that the in- 
ferred velocity dispersion of a galaxy is proportional to 
the square root of the mass of a galaxy over its size. 
At low redshift, where inferred velocity dispersion is cal- 
ibrated, stellar mass is a good estimate of the bary- 
onic mass content of a galaxy, especially for galaxies 
with high enough inferred dispersions to be included in 
this analysis. However, a t higher redshifts galaxies have 
highe r specific SFR s (e.g.lD addi et al.l l2007i: l Elbaz et al " 



2067 UNoeske et al.l l2007: Zh eng et al.ir 2007: Da men et al . 
120091) . If we assume the local IKennicutt, (,19981 ) relation 
between SFR and gas surface density holds, this implies 
that the gas density and therefore the gas mass fraction 
for a given size should increase at higher redshifts (e.g . 
IDaddi et al.ll2007l: iFranx et al.ll2008t INelson et al. l l2012l) . 
Assuming that the stellar mass of galaxies included in 
the calibration of inferred dispersion at z = will in- 
clude stars that have formed more recently, it could be 
important to include the gas mass in the total mass bud- 
get in estimating inferred velocity dispersions at earlier 
times. 

For each galaxy, we have estimates of the SFR from 
the SED fits. Using these values and assuming that gas 
is uniformly distributed in a circ ular aperture of r adius 
Te, we calculate gas masses using IKennicutt! (jl998D : 
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Fig. 1. — Rest-frame U-V and V-J colors selection diagrams for galaxies in the UDS and NMBS Cosmos fields for two redshift ranges. Grey 
scale boxes indicate the overall distribution of galaxies in the redshift range and colored points highlight galaxies above the cri„f > lOOkm/s 
threshold: red and blue points mark quiescent and star-forming galaxies. Median errors in rest-frame colors for selected galaxies are shown 
in the lower right corner of each panel. 

We calculate inferred dispersions using Mharyon = M^, + 
Mgas in place of Mi,. 

Although implied gas fractions increase for the cata- 
logs as a whole with redshift, gas fractions are relatively 
low, even for star-forming galaxies selected in this way. 
At lower redshifts, {Mgas /Mharyon) ^ 0.05 for galaxies 
above the velocity dispersion limit of the sample and 
the average does not evolve strongly with redshift. This 
is much lower than the 35% gas fractions measured by 
iTacconi et all (|20Tot ) in massive galaxies at z ^ 1.2. 
This discrepancy could be due to the fact that if ve- 
locity disper sion is a good pre dictor o f quiescence (e.g . 
IFranx et al.l 120081: IWake et al.l I2012hl : IBeh et all mm . 
then star-forming galaxies selected by velocity disper- 
sion would, on average, have lower SFRs than those with 
lower velocity dispersions, and therefore lower implied 
gas masses; an effect that would be magnified by the 
fact that the velocity dispersion cut increases with red- 
shift. This is consisten t with the fact that many of the 
z ^ 1.2 galaxies in the ITacconi et all (|2010f ) sample are 
consistent with being large, rotating disks with low ve- 
locity dispersions as opposed to those at higher redshifts, 
which might be more turbulent. Although we are not di- 
rectly measuring velocity dispersions, such massive star- 
forming galaxies have ver y large radii (rp ^ 5kpc) an d 
low Sersic indices (n < 1) |F6rster Schreiber et al.|[2011[ ). 



Even with M^ = 3 x 10^° Mq (|Tacconi et al.ll2010[ ). this 
corresponds to cri„/ ~ 80km/s, below the limits of this 
study. Additionally, this discrepancy could be due under- 
estimates of the star- formation rates calculated by stellar 
population synthesis modeling alone due to underestima- 
tion of dust attenuation, and the effect could be slightly 
stronger using SFRs based on MIPS fluxes. We note 
that some galaxies with high implied gas fractions are in 
our ffinf > lOOkm/s sample: in the highest redshift bin 
the most highly star-forming ^ 4% of (cin/ > lOOkm/s) 
galaxies in the Cosmos field and ~ 6% of galaxies in the 
UDS have average implied gas fractions equal to ~ 35%. 
We find that including the gas mass has very little effect 
on our final number density distributions (see §3) with 
no changes to the distribution of quiescent galaxies. At 
z > 0.6, there is a slight shift in the VDF for star-forming 
galaxies - at most a ^ 10% effect at a 200km/s, how- 
ever this shift is within the error bars. Given the uncer- 
tainties in the SFRs, we neglect the effects of cold gas 
on the inferred velocity dispersions for the rest of this 
analysis. For completeness, we include the VDF for star- 
forming galaxies with and without the inclusion of gas as 
Figure [TT] in Appendix D. 



2.3. Surface Densities 
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In addition to velocity dispersions, we calculate surface 
densities: 

^ = ^ (4) 

for all galaxies. This measure of galaxy compactness 
exhibits similar, although slightly weaker, corr elations 
with galaxy color (e.g. Fr anx et al, 20 08| : ^ Wake et al 
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and dark matter properties (jWake et al.ll2012al) 
as velocity dispersion. Therefore we will also express 
our results in terms of surface density in §3.3 to provide 
an important check on the results of this study. Inferred 
velocity dispersions, at least those > lOOkm/s, are well- 
calibrated for galaxies of all types in the SDSS. While 
there have been efforts to measure galaxy dynamics and 
determine the a ccuracy of i n ferred dispersion, at high 
redshifts Ce.g. iTreu et al.l l2005t Ivan der We i et al " 
20051 ICappellari et all I2009| : Ivan Dokkum et al 



2009t iNewman et al 



Bezanson et al. 2011 
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van de Sande et al.ll201lD these efforts have focused on 



quiescent galaxies with high inferred velocity disper- 
sions. We emphasize that we expect the results using 
surface density (oc M/B?) to be similar to those derived 
using inferred velocity dispersion (oc M/R), however, it 
is a more direct observational diagnostic and requires no 
calibration and is therefore not subject to uncertainties 
in evolution of constants unlike velocity dispersion. 

2.4. Differentiating between Star-Forming and 
Quiescent Galaxies 

We separate into red and dead versus blue and star- 
forming galaxies based on rest-frame colors. Although 
locally this would translate into an early type/late type 
split, this looser definition allows us to separate based 
on the stellar populations of galaxies while including the 
flexibility for structural evolution of entire populations 
of galaxies. 

Using EAzY, we calculate rest-frame U-V and V-J col- 
ors for each galaxy. Galaxies can have red U-V colors due 
to older stellar populations or because they are dusty. 
In the case where dust is obscuring the rest-frame UV 
and optical light, light is re-radiated at longer wave- 
lengths, causing these galaxies to have redder V-J col- 
ors than truly quiescent galaxies. We therefore create a 
two color plot with rest-frame U-V and V-J colors to iso- 
late red galaxies which are no longer forming stars from 
star-forming ga l axies, using the technique described by 
iWilliams et al.l (|2009|) . This color separation has been 
shown to agree well with other star forma tion indicators, 
such as MIPS 24ji t m detection (see e.g . IWilliams et al.l 
[2009; Wu vts et al.i r2009: Bell ct al. 2011) The thresholds 
drawn in the color-color plots to isolate the red sequence 
from star-forming galaxies evolve with time as the stellar 
populations of dead galaxies age. We adopt the re l ations 
determined for the NMBS fields in lWhitaker et al.l ()201lD 
in two redshift bins (0.3 < z < 0.5 and 0.5 < z < 1.5) and 
show the distributions in each field and redshift range in 
Figure [TJ with median random errors indicated in the 
lower right corner. The overall galaxy distribution in a 
given redshift range is indicated by the greyscale density, 
including galaxies below the velocity dispersion thresh- 
old, which we do not include in our analysis. Colored 
points indicate galaxies with ainf > lOOkm/s: star- 



forming galaxies are identified with blue points and qui- 
escent with red. In each plot, the red sequence of quies- 
cent galaxies is characterized by a diagonal sequence in 
the upper left region of color-color space. We note that 
although we define a bimodal population, there is a con- 
tinuum of galaxy ages and ongoing SFRs. Based on the 
correlation between velocity dispersion and galaxy type, 
we expect that even the star forming galaxies in the sam- 
ple may not be strongly forming new stars due to the cut 
in velocity dispersion, as suggested by the minimal im- 
pact of including gas mass (as discussed in §2.1). 

We evaluate the effect of defining quiescence using 
these color cuts by repeating the analysis using cuts on 
star-formation rates derived from IR and UV luminosi- 
ties in the Cosmos field, where MIPS 24 iim data is 
included in the NMBS catalogs (jWhitaker et al.l 120111: 
IWhitaker et al.l 120121) . We make a quiescent distinc- 
tion cut for star formation rates greater than twice the 
standard deviation below the redshift- dependent star- 
formin g main sequence as defined by IWhitaker et al.l 
()2012|) . Measured velocity dispersion functions for star- 
forming and quiescent galaxies using this method agree 
extremely well with those based on color cuts, after ap- 
plying a correction for the 24 /xm incompleteness as a 
function of a. The color derived quiescent population is 
less than 10% contaminated by MIPS detected sources, 
and down to 1 — 2% at the highest and lowest redshifts, 
corresponding to deviations from the observed VDF by 
an average of < 0.1 dex at all velocity dispersions. Be- 
cause of the additional complexity of understanding in- 
completeness in MIPS fluxes, we focus all further analysis 
on color-based quiescence as it is complete for all galax- 
ies in the surveys. We note that this test is performed 
in the NMBS Cosmos field, which may have more accu- 
rate rest-frame colors due to the improved sampling of 
the NIR medium band filters. There is a hint in Figure 
[T] that the two populations of galaxies are less distinct 
at higher redshifts in the UDS field than in the NMBS 
field and therefore the quiescent sample may be slightly 
more contaminated. We emphasize that although there 
are field-to-field offsets in the overall VDF, there does 
not appear to be a systematic trend in the relative frac- 
tions of star-forming and quiescent galaxies as a function 
of dispersion. 

In order to emphasize the importance of using velocity 
dispersion as opposed to stellar mass as an indicator of 
stellar properties, we plot the inferred dispersion versus 
steUar mass in Figure [Din the COSMOS field for 0.3 < 
z < 0.6 where the NEWFIRM medium band near-IR 
filters begin to contribute to the excellent photometric 
redshifts. It is clear from the histogram on the right 
that velocity dispersion produces a stronger separation 
between quiescent and star-forming galaxies than stellar 
mass (bottom histograms). 

3. NUMBER DENSITIES 

3.1. Velocity Dispersion Function for Star-Forming and 
Quiescent Galaxies 

Once we have calculated inferred dispersions and cat- 
egorized every galaxy as quiescent or star-forming, we 
calculate the number density as a function of velocity dis- 
persion for each field and population separately, included 
as open symbols (diamonds for the COSMOS field and 
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Fig. 2. — Inferred dispersion versus stellar mass in Figure [2] in 
the COSMOS field for 0.3 < z < 0.6, red points represent quiescent 
galaxies and blue points star-forming galaxies. Histograms on the 
bottom and right side indicate the relative contributions of the 
two populations. It is clear from these histograms that velocity 
dispersion of a galaxy is a better predictor of quiescence in galaxies 
than stellar mass. Approximate errors are indicated in the lower 
right corner. 

triangles for the UDS) in Figure [31 There is some field- 
to-field variation, however the qualitative distribution of 
galaxies is the same for the COSMOS and UDS fields. 
The volume-weighted average VDFs for quiescent and 
star-forming galaxies are shown in red and blue polygons. 
Error bars on all points reflect Poisson errors in addition 
to systematic errors due to cosmic variance, photometric 
redshifts and stellar masess (see §3.2). Incompleteness in 
Uinf is dominated by galaxies that are too faint to mea- 
sure reliable sizes and are therefore not included in the 
size catalogs [K > 22.4 in the UDS and K > 22.0 in the 
COSMOS held). We define the dispersion completeness 
limit (grey regions in Figure |3|) for each field and redshift 
range as the dispersion at which the 95% completeness 
plus 0.06 dex scatter about a linear fit to the K magni- 
tude - inferred dispersion relation reaches the magnitude 
limit of the size catalogs (see Appendix A). 

At every redshift, the shape of the quiescent VDF is 
very similar to the modified Schechter function that is 
fit to the V DF of elliptical galaxies in the SDSS (e.g . 
Sheth et al.l 1200 3: Mitchell et al.1 120051: IChoi et al.l 120071: 
Shankar et alBoOft .Bernardi et al.ll2010l ). On the other 
hand, the VDF for star-forming galaxies is primarily 
fiat in log-log space and is much steeper than the qui- 
escent VDF except at the highest dispersions. At all 
redshifts, star forming galaxies are much less common 
than quiescent galaxies at velocity dispersions greater 
than ~ 150km/s. 

In this analysis we construct the VDFs from inferred, 
not measured velocity dispersions. Therefore, the un- 
derlying distribution of intrinsic velocity dispersions is 
broadened by scatter in the ctq vs. amf relation, which 
becomes especially dominant at the steep, high disper- 
sion end of the VDF. This is essentially a form of Edding- 
ton bias whereby high ainf galaxies are likely to have 
intrinsically lower ctq values that have been scattered up- 
wards. We account for this by simulating an underlying 



distribution following a modified Schechter function for 
the quiescent galaxies and a double power law for the 
star- forming galaxies, scattering by the observed scatte r 
of 0.06 dex from the SDSS (see iBezanson et al.l[20Tll ). 
and then fitting the scattered distribution to the observed 
VDFs. The underlying distributions are included as blue 
(star-forming) and red (quiescent) dashed lines in Figure 
[31 This scatter should be interpreted as a minimal scat- 
ter, as additional sources of uncertainty such as redshift 
errors and systematic errors in stellar masses might in- 
crease with redshift. In every case, the unscattered VDF 
agrees with the scattered VDF at the low dispersion end. 
However, since the number density of very high disper- 
sion objects is low, adding scatter to the velocity disper- 
sion in the model broadens the intrinsic VDF. Given an 
increased scatter above the nominal 0.06 dex from SDSS 
the broadening of the intrinsic VDF would be stronger, 
implying an even steeper drop in the number density of 
galaxies with high velocity dispersions. This effect is rel- 
evant when interpreting the measured properties of small 
samples of high redshift galaxies selected on inferred ve- 
locity dispersion. 

To highlight the evolution of each population with 
time, we plot the VDF for each galaxy type at all red - 
shifts in the top row of Figure lil lBezanson et al.l (|2011[ ) 
demonstrated that the shape of the overall VDF has 
evolved strongly at the low velocity dispersion end, but 
that the overall number density of galaxies at high in- 
ferred velocity dispersions evolved very little since z ^ 
1.5. By separating into quiescent and star-forming galax- 
ies, we can gain insight into how this occurred. In gen- 
eral, the shape of the quiescent galaxy VDF remains sim- 
ilar in all redshift bins, but the low velocity dispersion 
end builds up with time. The number of star-forming 
galaxies at fixed dispersion is remarkably constant at all 
redshifts, perhaps with a slight increase at high disper- 
sions in the highest redshift bin. Another key feature is 
that there is very little evolution in the number of high 
velocity dispersion galaxies of either type. The strong 
evolution of quiescent galaxies and roughly stable num- 
ber density of star forming galaxi es has also been studie d 
as a function of stellar mass (e.g. iBrammer et al.ir2011[ ). 

3.2. Uncertainties in the Velocity Dispersion Function 

There are clear discrepancies between the number den- 
sities of galaxies in the two fields, often much larger 
than the associated Poisson errors. We will here try 
to evaluate the impact of these sources of error on the 
overall VDFs and include all error estimates along with 
the measured velocity dispersion functions in Table [1] 
We take the following sources of systematic error into 
account: cosmic v ariance estimated using QUICKCV 
(jMoster et al.l 12011V photometric redshift uncertainties 
and errors in stellar mass estimates. Because the UDS 
field relies only on broadband photometry, we also in- 
vestigate the possible influence of random errors on the 
photometric redshifts for the held. 

First, we we a dopt the metho dology of 
'S omerville et al.l (j2004D : IMoster et al.l (j2011[ ) to estimate 
the effects of cosmic variance on the me asured VDFs 
from individual fields. We use QUICKCV ([Moster et al.l 
1201 If ) to estimate the dark matter variance based on 
the volume probed by each redshift bin and survey and 
estimate the bias at each velocity dispersion based on 
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Fig. 3. — The VDF for quiescent (red) and star-forming (blue) galaxies as a function of redshift. Data points represent the measured 
VDF in Cosmos (diamonds) and UDS (triangles) fields individually. Large, light-colored regions show the combined VDF for the two 
galaxy populations at each redshift. Estimated incompleteness limits at a given redshift range for each photometric field are included as 
grey shaded regions. Here we see at all redshifts that the VDF for quiescent galaxies turns over at an intermediate velocity dispersion and 
falls off steeply at high inferred velocity dispersions. In contrast, the number density as a function of velocity dispersion is roughly flat for 
star- forming galaxies in all redshift bins. 



the b ias-number density relation in iSomerville et all 
(pool . In each case, the Poisson errors are larger than 
those expected by cosmic variance alone. However, 
we note that the difference between the two fields is 
often much greater than 1 — a cosmic variance and 
Poisson statistics would predict; the number densities 
of galaxies in the Cosmos field are higher than those in 
the UDS field in all but the highest redshift bin. The 
NMBS Cosmos field appears to be an overdense field, 
especially relative to the larger UDS field. We include 
the redshift distribution of galaxies in the surveys (above 
the minimum cri„/ > lOOkm/s threshold to be included 
in this study) in Figure 1101 In this figure, the large 
overdensities of galaxies at z ~ 0.72 and z ^ 0.88 in 
the Cosmos field are immediately apparent relative to 
the panel to the right showing the number density of 
galaxies in the UDS field in the same redshift range. 
These peaks fall in the same VDF redshift bin and can 
therefore explain the maximum of field to field variation 
observed at those redshifts. 

The two surveys included in this analysis are both ex- 
tremely deep, well-calibrated samples of galaxies, how- 



ever the NMBS Cosmos survey includes photometry 
in near-IR medium band filters that are explicitly de- 
signed to sample the Balmer/4000A break for galax- 
ies at z > 1, allowing for much better sampling of the 
SED and more accurate photometric redshifts. There- 
fore, while the UDS field should reach ~ 3% accuracy in 
^ (jWilhams et al.ll2009D . the NMBS field reaches an 

impressive ~ 1%— 2% accuracy. Photometric redshift un- 
certainties effect both the redshift bins in which galaxies 
fall, but also fold into their stellar masses and physical 
sizes. To determine the impact of systematic errors in 
redshifts on the final VDFs, we calculated the scatter 
in the VDF due to introducing redshift shifts to each of 
the fields, ±1% for Cosmos and ±3% in the UDS. We 
note that random errors in photometric redshifts would 
introduce a much more subtle effect as photometric red- 
shifts would scatter in either direction and would produce 
smaller error bars than this extreme estimate. The ef- 
fect of systematically shifting redshifts on the dispersion 
function is generally small, < 0.1 dex, comparable to the 
Poisson errors on each field. This is not the case, for ex- 
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Quiescent Galaxies Star-Forming Galaxies 




Fig. 4. — Top Row: VDFs for Quiescent (left panel) and Star- Forming (right panel) Galaxies at all redshifts. We note that while the 
quiescent VDF evolves strongly, building up low dispersion galaxies with time, the shape of the VDF for star-forming galaxies remains 
mostly constant. Bottom Row: Surface density functions for quiescent and star-forming galaxies with Poisson error bars. The differential 
evolution of number density for quiescent galaxies is also apparent as a function of surface density, as is the roughly constant shape for 
star-forming galaxies, at least out to z ~ 1.2. 



ample, in the Cosmos field between 0.6 < z < 1.2, where 
the errors grow larger (± ~ 0.2 — 0.3 or more). This 
can be explained by the tail of the z ~ 0.9 overdensity 
shifting between the two redshift bins. 

Systematic shifts in photometric redshifts will have a 
stronger impact on the overall VDF than random errors, 
however if the 3% accuracy of photometric redshifts is 
an underestimate in the UDS, these could become more 
important. To test this effect, we adopt a high estimate 
for random photo-z uncertainty, introducing a gaussian 
scatter the photometric redshifts with a width of 8% in 
Az/{1 + z), and refit stellar masses and calculate the 
final VDF. We verify that the effect of this test is quite 
minimal, introducing a mean offset of 0.04 dex (with an 
rms of 0.065 dex) for the quiescent VDF and 0.05 dex 
(rms=: 0.169 dex) in the VDF for star-forming galaxies. 
In every case the offset is less than the total uncertainty 
on the VDF for the field as quoted in Table 1. We con- 
clude that the inclusion of systematic errors in photo-zs 
provides a sufficient description for redshift uncertain- 
ties in the UDS. We note that this effect could indeed 
be increasing scatter in measured quantities and further 
broadening the high dispersion end of the VDF, even 



though we do not see evidence that the VDF in the UDS 
is shallower at the high velocity dispersion end than in 
the NMBS Cosmos field. 

Stellar masses are derived from SED fitting in both sur- 
veys, a method that is riddled with systematic sources 
of uncertainty. We adopt a 0.1 dex error in the stellar 
mass to allow for systematic uncertainties due to as- 
sumptions about IMF, star-formation histories and stel- 
lar population synthesis templates. We recalculate the 
velocity dispersions and velocity dispersion functions for 
logA/^ ± 0.1 dex. This has a very direct effect on the 
number densities of galaxies as they shift horizontally as 
galaxies shift to lower, and higher, bins in log(Ti„/, more 
strongly broadening the high dispersion end of the VDFs. 
These errors are nearly uniformly larger than other er- 
rors, especially for the VDF of star-forming galaxies, al- 
though the effect is slightly minimized near the peak of 
the quiescent VDF at (7i„f ^ 2.3, just due to the shape 
of the distribution. 

Altogether, the error estimates would add in quadra- 
ture and the individual fields are combined using a 
volume- weighted average. At most redshifts and velocity 
dispersions systematic errors in stellar mass dominate the 
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total error budget on the dispersion function. We include 
the total errors for every field as error bars on the indi- 
vidual points in Figure [3] Shaded regions in Figures [3l 
top panels of 21 and [8] reflect the volume- weighted VDFs 
with the combined total errors of the complete sample. 
By including these various sources of systematic error, 
the field to field variation of the VDF has become con- 
sistent with the total errorbars, which suggests that while 
we cannot discriminate amongst these sources of error, 
at least some of them must be important for proper in- 
terpretation of the VDF as measured with this data. 

3.3. Quiescent Fraction as a Function of Velocity 
Dispersion 

Another way of quantifying the build up of low disper- 
sion quiescent galaxies is to plot the quiescent fraction 
of galaxies as a function of inferred velocity dispersion, 
shown in the left panel of Figure [S] Unsurprisingly, the 
fraction of quiescent galaxies increases with velocity dis- 
persion in each redshift bin. However, an important re- 
sult of the buildup of low-dispersion quiescent galaxies 
and the constant star-forming VDF is that the disper- 
sion above which quiescent galaxies dominate becomes 
lower with time: from ~ 160km/s in the 1.2 < z < 1.5 
down to the lowest velocity dispersions included in this 
analysis of ^ 120 km/s by 0.3 < z < 0.6. This result sup- 
ports a picture of downsizing of galaxies, with galaxies 
turning off at lower dispersions, i.e. larger radius and/or 
smaller mass, and building up the red sequence at later 
times. We note that this trend is most significant when 
the systematic errors in stellar mass cancel out in calcu- 
lating quiescent fractions; black error bars reflect Pois- 
son, cosmic variance and systematic redshift errors only, 
grey error bars in Figure [Sb include all errors discussed in 
§3.2. This result is consistent with the observational re- 
sult that quiescent galaxies with low disper sions have sys- 
tema t ically younger ages i n the SDSS (e.g. lHaiman et all 
[2007t iShankar et all [20091 ). The right panel of Figure \5\ 
shows the threshold velocity dispersion as a function of 
redshift and a fit such that athresh oc (1 -I- z)^". 

3.4. Surface Density Function for Star-Forming and 
Quiescent Galaxies 

In the local Universe, inferred velocity dispersion is 
well determined by stellar mass and structural param- 
eters of galaxies - regardless of morphology or star- 
formation rates. However, all measurements of the veloc- 
ity dispersions of galaxies at high redshifts are made for 
quiescent galaxies and inferred velocity dispersion is es- 
sentially uncalibrated for star-forming galaxies at z > 0. 
We therefore investigate whether the number density as 
a function of surface density, calculated in §2.4, which 
may be a more observationally motivated quantity for 
star-forming galaxies, exhibits similar evolutionary be- 
havior as the inferred VDF. We present the "surface 
de nsity function" (see also the "compactness function" 
in I Newman et al.] ([2012)) for quiescent and star-forming 
galaxies at all redshift bins in the bottom row of Figure 
m emphasizing that the overall shapes and evolution of 
these distributions are quite similar to the VDFs (shown 
in the top row). Symbols and colors are the same as 
for the VDFs and completeness limits are derived as de- 
scribed in §2.1 by using the scatter about the K magni- 
tude - surface density relation. 



We note that while the surface density function is an 
important test of these results, the Sersic dependence of 
inferred velocity dispersion could make the VDF more 
relevant for comparison with simulations because it can 
account for the structural evolution that has occurred 
in ga l axies in this redshift range (e.g. Ivan der Wei et ahl 
120111: iBuitrago et all [20110 . Furthermore, while there 
may not be physical constraints on the rate at which 
galaxy surface densities can evolve with time, the mech- 
anisms that can increase or decrease the motion of stars 
in a galaxy are more constrained. 

4. A SIMPLE EXPLANATION FOR THE COEVOLUTION 
OF THE QUIESCENT AND STAR-FORMING VDFS 

4.1. Model Overview 

As we discuss in §1, the velocity dispersion of an indi- 
vidual galaxy can also evolve and we can make some sim- 
plistic statements about the physical processes that could 
be responsible for these changes. Gas dissipation would 
increase the velocity dispersion of a galaxy, which would 
be primarily important for star-forming galaxies. Addi- 
tionally, the contribution of dark matter to the dynam- 
ical mass of a galaxy could evolve with redshift, which 
would alter the ratio between the stellar and dynami- 
cal mass for a galaxy. This would therefore change the 
galaxy's inferred velocity dispersion just due to evolution 
of the dark matter halo. Merging of galaxies could either 
increase or decrease the velocity dispersion of a galaxy, 
depending on the initial orbital energies of the progenitor 
galaxies and their relative masses and gas content. For 
example, it has been suggested that minor merging could 
puff up galaxies, thereby decreasing the mass within 
a given radius and decreasing velocity dispersion (e.g . 
iBovlan-Kolchin et al.' '2005', '2006"; 'Bezan son et"all 120091 : 
^aab et al. 2009; Oscr ct al. 2011). Finally, gas outflows 
could decrease the velocity dispersion of a galaxy as mass 
is expelled from the central regions of the galaxy. In ad- 
dition to changes in the velocity dispersion of individual 
galaxies in these samples with time, the number density 
of the star-forming and quiescent galaxies can change due 
to transformation between the populations. 

In order to demonstrate the qualitative effects of these 
processes on the quiescent and star- forming VDFs, we 
produce a simple model in which the number density 
of galaxies initially follows the underlying VDFs fit to 
the 1.2 < z < 1.5 redsh i ft bin and evol ves for ward in 
time (also s ee iBell et al.l (|2007[ ): Walclie r et al.l (2008); 
IPeng et al.l (pOlfli r for similar analysis based on the evo- 
lution of the stellar mass function). Without assign- 
ing physical mechanisms as responsible for changes in 
the galaxy population, we allow individual star-forming 
galaxies to quench and transform into quiescent. While 
this simple evolution could reproduce the build-up of the 
quiescent VDF, it would deplete the population of star- 
forming galaxies and contradict the observed constancy 
of the VDF for star-forming galaxies. Therefore, in order 
to maintain a constant distribution of star-forming galax- 
ies, we allow individual galaxies to undergo net shifts in 
velocity dispersion in a redshift interval, d\oga/dt. The 
effects of these two processes on the VDFs are indicated 
in Figure [51 In this case, transformation moves galax- 
ies from the blue star-forming VDF to the red quiescent 
VDF, decreasing the number density of the former and 
increasing the latter. Net increases in velocity dispersion 
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Fig. 5. — Left Panel: Fraction of quiescent galaxies as a function of inferred velocity dispersion. The dashed black line indicates the 
50% threshold between quiescent or star-forming dominance. At all redshifts, the quiescent galaxies dominate at high velocity dispersion, 
however the crossover at which quiescent galaxies become more prevalent than star forming galaxies evolves to higher velocity dispersions 
at higher redshift. Right Panel: Evolution of the threshold velocity dispersion between quiescent and star-forming galaxies as a function 
of redshift. This cross-over velocity dispersion decreases with time to the power (Jthresh oc {1 + z)^", as denoted by the red line. 

ative. Adding only transformation to the model quickly 
depletes the number of star-forming galaxies, whereas 
the observed star-forming VDF is constant with time. 
Therefore, we introduce a net change in the velocity dis- 
persions of star-forming galaxies at a given velocity dis- 
persion, such that the number of galaxies leaving the bin 
in velocity dispersion equals the number of galaxies en- 
tering from the adjacent bin. Matching the number flux 
insures the constant shape of the star-forming VDF with 
time. At each timestep, the velocity dispersions of in- 
dividual galaxies are randomly scattered by 0.06 dex to 
reflect the scatter in measured and inferred velocity dis- 
persions in the SDSS and the VDF is recalculated. 

Figure [7] demonstrates the overall properties required 
for the model to reproduce the observed VDF evolu- 
tion. In the left panel of Figure [7] we present the im- 
plied growth fraction of quiescent galaxies as a function 
of velocity dispersion between each pair of redshift bins. 
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Fig. 6. — Simulated VDFs from the highest redshift bin with ar- 
rows illustrating the effects of transformation and dissipative pro- 
cesses on the shape of the VDF. 

would shift the star- forming VDF to the right. 

4.2. Model Fitting 

In this context, the increase in low-dispersion quiescent 
galaxies can simply be explained by the transformation 
of star-forming galaxies where they are most common - 
at the low dispersion end of the VDF. We set the trans- 
formation rate equal to the increase in the observed VDF 
for quiescent galaxies at a given velocity dispersion and 
in each redshift bin as determined by the unscattered fits. 
We impose a minimum transformation rate of zero for ve- 
locity dispersions at which the implied evolution is neg- 



f .\- Nq{cr,z)-Nq{(7, Zj) 
Jgrowth(Cr,Z^ — [0 

Nq{cr,z) 

As suggested by the evolution of the observed VDF, the 
growth in the number density of quiescent galaxies is 
strongest at lower inferred dispersions. The opposite 
trend is exhibited by the quenching fraction (top right 
panel) of Figure U\ 



f quench (^^7 ^) 



Ngja, z) ~ Ngja, Zj) 
Nsf{c7,z) 



(6) 



At the low velocity dispersions, star-forming galaxies 
are most common and the quenched fraction is low. 
At high velocity dispersions, a > 300km/s, the VDF 
does not evolve strongly and therefore the galaxies re- 
main unchanged. At intermediate velocity dispersions 
(cr ~ 250km/s), where few galaxies are quenching in 
a given timestep, the reservoir of star-forming galaxies 
is also minimal, therefore the quenched fraction is the 
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Fig. 7. — Top Left Panel: Fraction of newly quenched galaxies relative to the existing number of quiescent galaxies as a function of 
velocity dispersion. At all redshifts, this fraction increases from zero evolution at the highest redshifts to ~ O.lperlOOMyr at low velocity 
dispersions. Top Right Panel: Fraction of star-forming galaxies that are transformed into quiescent galaxies as a function of velocity 
dispersion. This fraction exhibits the opposite behavior, with intermediate dispersion galaxies quenching much more rapidly than those 
with lower dispersions. Again there is no evolution for galaxies with the highest dispersions. Bottom Left Panel: Rate of change in the 
velocity dispersions of star-forming galaxies as a function of velocity dispersion necessary to sustain the observed galaxy transformation and 
maintain a constant star-forming VDF. This model requires extremely rapid changes in the velocity dispersions of star-forming galaxies, 
increasing with velocity dispersion. Bottom Right Panel: Median redshift at which galaxies quench in the simulation as a function of 
velocity dispersion. In the model, high dispersion galaxies quench first and galaxies with lower velocity dispersions quench at later times. 



highest at these velocity dispersions. Finally, the rate at 
which star-forming galaxies enter into a bin of velocity 
dispersion implies an average rate of change in velocity 
dispersion, which is shown as a function of inferred dis- 
persion in the bottom right panel of Figure [T] 

/dlogcr \ fquench{T, Z) 

\^("'"V = ^""^ ^ — It — 

Here it is clear that the increasing quenching fraction 
with velocity dispersion implies that the rate at which 
the dispersions of galaxies must grow to maintain a con- 
stant distribution also increases with velocity dispersion. 
This all combines to create a picture of transformations 
of intermediate dispersion star-forming galaxies, imply- 
ing that at least in this regime, the quenching mechanism 



is rapid. At the same time, some process (or set of pro- 
cesses) efficiently drives mass to or creates mass in the 
centers of galaxies more rapidly than growing their sizes, 
thus increasing their inferred velocity dispersions and re- 
plenishing the ranks of star-forming galaxies at a given 
velocity dispersion. 

4.3. Model Results and the "Reservoir Problem" 

The results of this simulation are included in Figure [51 
which demonstrates that the evolution of the quiescent 
and star-forming VDFs can be reproduced by this illus- 
trative model. The high velocity dispersion end of both 
VDFs is roughly constant in the simulation and there 
is a gradual build up of quiescent galaxies with low ve- 
locity dispersions. As a result, quiescent galaxies with 
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Fig. 8. — Evolution of the observed (large, transparent polygons) quiescent {Left Panel) and star-forming [Right Panel) VDFs and those 
from simple model (narrow, opaque, outlined polygons) including transformation and evolution in the velocity dispersions of individual 
star-forming galaxies. Transformation rates are determined from the growth of the observed quiescent VDF. In order to maintain the 
constant VDF for star- forming galaxies, dloga/dt for each bin in velocity dispersion is tuned such that the number of galaxies increasing 
dispersions out of the bin matches the number of quenched galaxies. 

galaxies as estimated by photometric stellar masses and 
sizes. We find three striking results: 



lower velocity dispersions will have the youngest stellar 
ages, which is consistent with obse rvations of galaxies in 
the SDSS (e.g. lHaiman et al1l2007| ) and paints a comple- 
mentary picture of the simulated b uild-up of the VPF fo r 
early type galaxies as presented in'S hankar et al.l ()2009[ ). 

It is clear that the high quiescent fraction and increas- 
ing number density of galaxies with intermediate veloc- 
ity dispersions requires rapid quenching and structural 
evolution of star-forming galaxies. While this result is 
quite extreme, it is hard to escape given the current 
observations. One potential caveat is that there could 
be some selection bias such that current surveys of high 
redshift galaxies systematically missing obscured, star- 
forming and dynamically massive galaxies or galaxy cen- 
ters which later evolve into compact quiescent galaxies. 

We emphasize that this model by no means encom- 
passes all physical processes that would in fact affect the 
evolution of the velocity dispersion functions. Additional 
complications could include evolution in the dispersions 
of quiescent galaxies, perhaps due to the effects of minor 
merging, which has been shown to be important in de- 
creasing the velo city dispersions o f galaxies in cosmolog- 
ical simulations (jOser et al.l 1201 If) . Furthermore, major 
mergers could affect the velocity dispersion of quiescent 
galaxies, depending on the initial relative galaxy orbits 
and equal mass merging would decrease number densi- 
ties. However, while these effects might shuffle the dis- 
tribution of quiescent galaxies, the qualitative agreement 
of our model should still apply. Even with these possible 
complications, it appears that the VDFs require a strong 
build up of quiescent galaxies at low velocity dispersions 
combined with rapid quenching and increasing velocity 
dispersion for star-forming galaxies with intermediate ve- 
locity dispersions. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we investigate the evolution of the Veloc- 
ity Dispersion Functions for star-forming and quiescent 



• At the highest velocity dispersions (> 300 km/s), 
the number density of galaxies evolves very little 
with time if at all for both quiescent and star- 
forming galaxies. This implies that galaxies with 
high velocity dispersions, like the halos they oc- 
cupy, form at early times and, while they may 
evolve structurally, their potential wells and inter- 
nal dynamics remain roughly unchanged with time. 

• The number density of quiescent galaxies with 
moderate to low velocity dispersions increases with 
time, presumably as star-forming galaxies with 
those dispersions shut off their star formation. 

• The VDF for star-forming galaxies remains roughly 
constant over the majority of cosmic time from 
0.3 < z < 1.5. While there is potentially concern 
about the accuracy of inferred velocity dispersions 
of star-forming galaxies at z > 0, these results also 
hold for the number density of star-forming galax- 
ies as a function of surface density out to z ^ 1.2. 

We proposed a simple model to reproduce these ef- 
fects including transformation of star-forming to qui- 
escent galaxies and increasing velocity dispersions for 
star-forming galaxies. In this case, galaxies with a > 
100 km/s must quench rapidly, with an efhciency that 
increases with velocity dispersion up to cr 250 km/s 
and rapid increases in velocity dispersions to maintain 
t he distribution of sta r-forming galaxies. 

IWake et al.l (|2012a[ ) demonstrated that the velocity 
dispersion, and to a slightly lesser extent the surface den- 
sity, of a galaxy is a good predictor of clustering and 
therefore possibly dark halo mass. Therefore, if this re- 
sult holds within the velocity dispersion range probed 
by this study and velocity dispersion remains a tied to 
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halo mass out to z ~ 1.5, we can speculate about the 
distribution of star-forming and quiescent galaxies as a 
function of velocity dispersion in the context of occupa- 
tion of dark matter halos. For example, given the stabil- 
ity of the distribution of velocity dispersions, or surface 
densities, for star-forming galaxies, these results would 
imply that the processes of accretion and gas cooling 
are tied to halo mass and at least in the most massive 
halos do not evolve with redshift, at least to z ~ 1.2. 
Given the evolution of the quiescent fraction, this would 
also suggest that star formation primarily turns off in 
lower mass halos at later times. However this simple 
picture is complicated by the fact that the velocity dis- 
persion (or size and stellar mass) of a galaxy can likely 
evolve more rapidly than the dark matter halo implies 
that this relation is likely to evolve with time, or at least 
exhibit a fair amount of scatter. Furthermore, the fact 
that star-formatio n rates at a given stellar mass increase 
with redshift (e.g.lDadd i et al.ll2007l: INoeske etllll2007l: 
IWhitaker et "81112011 " contradicts a picture of constant 
gas accretion and cooling with time. 

Overwhelming evidence has shown that every galaxy 
hosts a SMBH in its center and that the mass of 
that black hole is well correlated w ith the velocity 
dispersio n of the host's bulge (e.g. iMagorrian et al.l 
1998; Fe rrarese fc MerrittI [2000: iGebhardt et al.l 120001: 
Giiltekin et al.l I2Q09 D. Assuming that this correlation 
does not evolve strongly with time, the redshift evolution 
of the VDF implies that the build-up of SMBHs is anti- 
hierarchical; the most massive SMBHs are in place by 
z ~ 1.5, whereas the smaller black holes are still growing 
at lower redshifts. This "down-sizing" is both consistent 
with other studies of the evolution of the SMBH function 
and with the idea that the most massive black holes are 
relics of the most active AGN activity, which peaked at 
z ~ 2. 

This work relies heavily on the accuracy of inferred 
velocity dispersions, which are calibrated locally. Our 
conclusions are strengthened by the fact that the num- 
ber density of galaxies as a function of surface den- 



sity exhibits quite similar behavior. However, all of 
this analysis relies on the accuracy of Sersic profile 
fits and stellar mass estimates. Size measurements are 
largely made fr om ground-based imaging and although 
iBezanson e t al. (2011) demonstrated the consistency of 
sizes in the NMBS COSMOS catalog with those mea- 
sured from WFC3 imaging, the results would be less un- 
certain with the use of exclusively space-based imaging 
in the rest-frame optical, especially for the most com- 
pact galaxies. Although corrections for non-homology 
via the Ky{n) are small, repeating this study using only 
space-based imaging would vastly improve our ability ac- 
curately estimate best-fit Sersic parameters. Addition- 
ally, size measurements made with deeper imaging would 
allow us to probe to lower inferred dispersions at the 
higher redshifts probed by this study. Stellar masses are 
based on a number of assumptions about the star forma- 
tion h istories of galaxies and IMF; we adopt a iChabrieil 
(|2003f ) IMF, however there have been hints that this 
could un der-predict masses for the highest mass galax- 
ies (e .g. Ivan Dokkum &: Conrovl [20ToriCappellari et al.l 
I20TI . Finally, the accuracy of inferred velocity disper- 
sion as a predictor of the actual dynamics of galaxies at 
z > remains an open question, especially for galaxies 
which are star-forming or with low velocity dispersions. 
This highlights the need for a more complete, statistical 
sample of measured velocity dispersions at z > to de- 
termine the accuracy of inferred velocity dispersions and 
evaluate the direct measurements of galaxy dynamics at 
z>0. 
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APPENDIX 

A. INCOMPLETENESS LIMITS IN VELOCITY DISPERSION 

Both the NMBS Cosmos and UKIDSS UDS photometric catalogs probe to remarkably faint magnitudes and therefore 
are complete to low stellar masses at all redshifts probed by this study (logAf^, ^ 9 at the lowest redshifts and 
log Mi, ^ 10 at z ^ 1.5, however size measurements were only made for a bright subsamplc of galaxies in each survey. 
For the NMBS Cosmos field, measurements down to a K magnitude limit of 22, for the UDS field K > 22.4. Therefore, 
incompleteness in this study will be dominated by galaxies that are too faint to be included in the size catalogs. In 
Figure[9]we show the K magnitude vs. inferred velocity dispersion for separately for each field and each redshift range. 
In each case, galaxies fall around a generic trend (as indicated by diagonal lines: red-UDS, green-Cosmos) with a fair 
amount of scatter, which depends on velocity dispersion. We fit the trend in each field for loga > 2.0. While the 
trend evolves to brighter magnitudes at low redshift, the relations measured agree extremely well from field to field. 
Assuming that these relations hold below the magnitude limits of our survey, we define conservative incompleteness 
limits as the velocity dispersion by calculating a 2 — cr upper limit to the K magnitude for a given velocity dispersion 
and calculate the limiting cr, plus 0.06 dex for the scatter between ctq and (Tm/, at the magnitude limit for each redshift 
and field (indicated by vertical lines in Figure |9|). Therefore, we shade regions in Figure [9] that are excluded from our 
analysis either due to the general logcr > 2.0 cut or our defined incompleteness limits. 

B. REDSHIFT DISTRIBUTION OF GALAXIES IN THE NMBS COSMOS AND UKIDDS UDS FIELDS 

Simple estimates of the cosmic variance for these redshift bins predict much smaller discrepancies than are observed 
between the UDS and Cosmos surveys. In Figure [10] we present the redshift distributions of galaxies in these two 
surveys. In particular we note the peaks in the Cosmos field, primarily the overdensity at z ~ 0.7 and at z ~ 0.9, 
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Fig. 9. — K magnitude vs. inferred velocity dispersion for UKIDSS UDS and NMBS Cosmos fields. Black points represent galaxies in 
a given redshift range, horizontal lines reflect the magnitude limits of the size catalogs (Red for UDS at K = 22, Green for Cosmos at 
K = 22.4). At each redshift incompleteness defined as the velocity dispersion at which twice the scatter about the overall linear trend 
(diagonal lines) meets the magnitude limit for the survey (horizontal lines). Regions excluded from our measured VDFs, either due to 
incompleteness or the cri„f > lOOkm/s limit, are shaded in grey. Median error bars are included in the lower right corner of each panel. 

which likely boost the number density of galaxies in the 0.6 < z < 0.9 redshift bin. 

C. VELOCITY DISPERSION FUNCTIONS AND ERROR ESTIMATES 

We include the measured velocity dispersion functions in tabular form with all calculated error estimates: Poisson 
errors (±Poi), cosmic variance (±ct;), systematic redshift uncertainties (±(z)), systematic uncertainties in stellar 
mass measurements (±(m)) and the total error bars with all errors added in quadrature for the individual fields and 
volume-weighted combined measurements for quiescent galaxies ($q) and star- forming galaxies ($s/)- 
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Fig. 10. — Redshift distributions of galaxies in the NMBS Cosmos and UKIDSS UDS fields. Black histograms reflect the overall 
distribution of galaxies and blue/red histograms show the distribution of star- forming and quiescent galaxies respectively. 

TABLE 1 — Continued 



logai„/ Field log*, ±{Poi)±(cv) ±(z) ±(m) ±{Tot) log<J>^j ±(Poi)±(cv) ±(z) ±{m) ±{Tot) 
[km/s] [Mpc~^dex"^] [Mpc~^dex~^] 







Total 


"2.55 


0.03 


0.04 


+0.05 
-0.01 


+0.14 
-0.30 


+0.16 
-0.30 




-3.55 


0.09 


0.05 


+0.01 
-0.16 


+0.08 
-0.12 


+0.14 
-0.23 


2.4- 


2.5 


































Cosmos 


-2.89 


0.09 


0.06 


+0.05 
-0.00 


+0.31 
-0.17 


+0.33 
-0.20 




-3.99 


0.31 


0.09 


+0.48 
-0.30 


+0.60 
-0.30 


+0.83 
-0.53 






UDS 


-3.34 


0.08 


0.06 


+0.07 
-0.11 


+0.33 
-0.48 


+0.36 
-0.50 




-4.77 


0.43 


0.11 


+0.30 
-99.00 


+0.85 
-99.00 


+ 1.00 
-99.00 






Total 


-3.18 


0.06 


0.05 


+0.06 
-0.08 


+0.26 
-0.37 


+0.28 
-0.39 




-4.42 


0.25 


0.07 


+0.26 
-99.00 


+0.56 
-99.00 


+0.69 
-99.00 


2.5 - 


2.6 


































Cosmos 


-3.33 


0.14 


0.07 


+0.05 
-0.65 


+0.28 
-0.65 


+0.32 
-0.94 


< 


-4.29 
















UDS 


-4.47 


0.31 


0.09 


+0.00 +0.65 
-0.00 -99.00 


+0.73 
-99.00 


< 


-4.77 
















Total 


-3.85 


0.13 


0.06 


+0.03 +0.31 
-0.39 -99.00 


+0.35 
-99.00 


< 


-4.89 












2.6 - 


2.7 


Cosmos 


< -4.29 












< 


-4.29 













Evolution of Quiescent and Star- Forming Galaxies Since 2; ~ 1.5 
TABLE 1 — Continued 



log ainf 

[km/s] 


Field 


log*q ±{Poi)±{cv) 
[Mpc-^dex-i] 


±(z) ±{m)±{Tot) loi 

[Mpc" 


l^sf ±iPoi)±icv) ±{z) ±{m) ±{Tot) 




UDS 


< -4.77 


< 


-4.77 




Total 


< -4.89 


< 


-4.89 


0.6 < z < 0.9 



2.0-2.1 



2.1 



2.2 



2.3 



2.4 



2.5 



2.6 



2.2 



2.3 



2.4 



2.5 



2.6 



2.7 



0.9 < z < 1.2 



2.0-2.1 



CoSTTiOS 


—2 


82 





06 





05 


+0.06 
-0.00 


+0.13 
-0.06 


+0.16 
-0.10 




-2 


51 





04 





05 


+0 
-0 


15 

19 


+0 
-0 


09 

10 


+0 
-0 


19 

23 


UDS 


-2 


79 





03 





04 


+0.00 


+0.06 


+0.08 




-2 


71 





03 





04 


+0 


05 


+0 


11 


+0 


13 


Total 


-2 


79 





03 





04 


+0.01 

Q 


+0.06 
Q Qg 


+0.07 
Q 




-2 


65 





02 





03 


+0 

Q 


05 
03 


+0 

Q 


10 
09 


+0 

Q 


12 
10 


(In ^Tfin^ 


—2 


62 





05 





05 


+0.09 
-0.03 


+0.12 
-0.07 


+0.16 
-0.10 




-2 


72 





05 





05 


+0 
-0 


14 
09 


+0 
-0 


12 
08 


+0 
-0 


19 
14 


UDS 


-2 


71 





03 





04 


+0.06 


+0.02 


+0.08 




-2 


93 





04 





04 


+0 


05 


+0 


13 


+0 


15 


Total 


-2 


69 





02 





04 


+0.06 


+0.02 


+0.07 




-2 


87 





03 





04 


+0 


05 


+0 


11 


+0 


13 


Cosmos 


—2 


35 





03 





05 


+0.03 
-0.09 


+0.05 
-0.15 


+0.08 
—0.18 




-2 


91 





06 





05 


+0 
-0 


04 
08 


+0 
-0 


11 

09 


+0 
-0 


14 

15 


UDS 


-2 


60 





03 





04 


+0.00 


+0.00 


+0.05 




-3 


39 





06 





05 


+0 


09 


+0 


24 


+0 


27 


Total 


-2 


52 





02 





03 


+0.00 


+0.01 


+0.04 




-3 


21 





04 





04 


+0 


07 


+0 


18 


+0 


20 


CoSTflOS 


—2 


31 





03 





05 


+0.08 
-0.11 


+0.02 
-0.05 


+0.10 
-0.13 




-3 


18 





09 





06 


+0 
-0 


08 
06 


+0 
-0 


18 
40 


+0 
-0 


22 
42 


UDS 


-2 


79 





03 





04 


+0.00 
—0.02 


+0.15 
—0.16 


+0.16 

—0.17 




-4 


06 





14 





07 


+0 
— 


15 

00 


+0 
— 


36 

15 


+0 

— 


42 

22 


Total 


-2 


61 





02 





03 


+0.02 
-0.03 


+0.11 
-0.12 


+0.12 
-0.13 




-3 


64 





07 





05 


+0 

-0 


09 
03 


+0 

-0 


21 
20 


+ 
-0 


25 
22 


Cosmos 


-2 


49 





04 





05 


+0.13 
-0.22 


+0.13 
-0.22 


+0.19 
-0.32 




-3 


80 





18 





08 


+0 

-0 


12 
30 


+0 

-0 


22 
30 


+0 

-0 


32 
47 


UDS 


-3 


32 





06 





05 


+ 0.05 
-0.01 


+ 0.37 
-0.59 


+ 0.38 
-0.60 




-4 


58 





25 





08 


+ 
-0 


00 
48 


+ 
-0 


37 
48 


+ 
-0 


45 
72 


Total 


-2 


94 





03 





04 


+0.06 
-0.09 


+0.22 
-0.36 


+0.24 
-0.37 




-4 


23 





14 





06 


+0 
-0 


05 
31 


+0 
-0 


24 
31 


+0 
-0 


30 
47 


Cosmos 


-3 


13 





08 





06 


+0.23 
-0.45 


+0.42 
-0.67 


+0.49 
-0.81 


< 




4.58 






















UDS 


-4 


76 





31 





09 


+0.30 
-0.00 


+0.85 
-99.00 


+0.95 
-99.00 


< 




5.06 






















Total 


-3 


71 





08 





05 


+0.20 
-0.37 


+0.37 
-99.00 


+0.44 
-99.00 


< 




5.19 






















Cosmos 


-4 


58 





43 





11 


+0.60 
-99.00 


+0.78 
-0.00 


+ 1.08 
-99.00 


< 




4.58 






















UDS 


< - 


5.06 
















< 




5.06 






















Total 


-5 


19 





43 





11 


+0.60 
-99.00 


+0.78 
-0.00 


+ 1.08 
-99.00 


< 




5.19 






















.2 


Cosmos 


-3 


21 





08 





06 


+ 0.19 
-0.31 


+0.14 
-0.10 


+0.26 

-0.34 




-2 


60 





04 





05 


+0 

-0 


08 

12 


+0 

-0 


07 
06 


+0 

-0 


12 

15 


UDS 


-3 


21 





04 





04 


+0.07 
-0.01 


+0.06 
-0.09 


+0.11 
-0.11 




-2 


74 





03 





04 


+0 
-0 


04 
00 


+0 
-0 


11 
13 


+0 
-0 


13 

14 


Total 


-3 


21 





04 





04 


+0.07 
-0.03 


+0.06 
-0.08 


+0.11 
-0.11 




-2 


70 





02 





03 


+0 
-0 


04 
02 


+0 
-0 


10 

11 


+0 
-0 


11 

12 


Cosmos 


-2 


92 





05 





05 


+0.17 
-0.18 


+0.11 
-0.14 


+0.22 
-0.24 




-2 


80 





05 





05 


+0 
-0 


01 
00 


+0 
-0 


15 
13 


+0 
-0 


16 
15 


UDS 


-3 


03 





04 





04 


+0.08 
-0.00 


+0.09 
-0.12 


+0.13 
-0.13 




-2 


96 





03 





04 


+0 
-0 


07 
00 


+0 
-0 


09 
13 


+0 
-0 


13 
14 


Total 


-3 


00 





03 





04 


+0.07 
-0.02 


+0.08 
-0.10 


+0.12 
-0.12 




-2 


92 





03 





03 


+0 
-0 


06 
00 


+0 
-0 


08 
11 


+0 
-0 


11 
12 


Cosmos 


-2 


81 





05 





05 


+0.16 
-0.17 


+0.00 
-0.01 


+0.17 
-0.18 




-3 


01 





06 





05 


+0 
-0 


05 
12 


+0 
-0 


08 
16 


+0 
-0 


12 
21 


UDS 


-2 


94 





03 





04 


+0.04 
-0.00 


+0.00 
-0.00 


+0.07 
-0.05 




-3 


30 





05 





05 


+0 

-0 


05 
00 


+0 

-0 


21 

31 


+0 

-0 


22 

31 


Total 


-2 


91 





03 





03 


+0.04 
-0.02 


+ 0.00 
-0.00 


+0.06 
-0.05 




-3 


21 





04 





04 


+0 
-0 


04 
02 


+0 

-0 


17 

25 


+0 

-0 


19 

26 


Cosmos 


-2 


73 





04 





05 


+0.11 
-0.12 


+ 0.00 
-0.09 


+0.13 
-0.17 




-3 


37 





09 





06 


+0 
-0 


09 
00 


+0 

-0 


21 

13 


+0 
-0 


25 
17 


UDS 


-2 


98 





03 





04 


+0.06 
-0.01 


+0.04 
-0.16 


+0.09 
-0.17 




-3 


98 





11 





06 


+0 

-0 


30 
00 


+0 

-0 


37 
33 


+0 

-0 


49 

35 


Total 


-2 


90 





03 





03 


+0.05 
-0.02 


+0.03 
-0.14 


+0.08 
-0.15 




-3 


73 





07 





05 


+0 
-0 


21 
00 


+0 
-0 


26 
23 


+0 
-0 


35 
25 


Cosmos 


-2 


90 





05 





05 


+0.04 
-0.09 


+0.14 
-0.15 


+0.16 
-0.19 




-3 


82 





15 





08 


+0 
-0 


00 
43 


+0 
-0 


33 
90 


+0 

-1 


37 
01 



2.1 



2.2 ■ 



2.3- 



2.4- 



2.2 



2.3 



2.4 



2.5 



16 



TABLE 1 



Continued 



log a, 


inf 


Field 


log*. 


±{Poi 


)±{cv) ±{z) ±{m)±{Tot) 


log $5 J 


±{Poi)±{cv) ±{z) ±{m)±{Tot) 


[km. 


/s] 


1 


'Mpc~^dex~^ 


1 










1 


'Mpc~^dex~ 










UDS 


—3.41 


0. 


05 


0.05 


+0.10 
-0.00 


+0.27 
-0.35 


+0.29 
-0.36 


—4.91 


0.31 


n no +0.54 +0.60 +0.87 
^•'-'^ -0.00 -0.30 -0.44 






Total 


—3.22 


0. 


,04 


0.04 


+0.07 
-0.02 


+0.20 
-0.26 


+0.22 
-0.27 


—4.33 


0.14 


n nfi +0.23 +0.32 +0.43 
U.UD _o.24 -0.53 -0.61 


2.5 - 


2.6 




























Cosmos 


-3.43 


0, 


,10 


0.06 


+0.10 
-0.00 


+0.37 
-0.22 


+0.40 
-0.25 


< -4.73 










UDS 


-4.09 


0, 


,12 


0.06 


+0.00 
-0.16 


+0.33 
-0.51 


+0.36 
-0.55 


-5.21 


0.43 


r, -I n +0.00 +0.00 +0.45 
U-IU _9g 00 -99.00 -99.00 






Total 


-3.81 


0, 


,08 


0.05 


+0.03 
-0.11 


+0.25 
-0.35 


+0.27 
-0.38 


-5.33 


0.43 


n ip, +0.00 +0.00 +0.45 
U-IU _gg 00 -99.00 -99.00 


2.6 - 


2.7 




























Cosmos 


-3.77 


0, 


,14 


0.07 


+0.00 
-0.65 


+0.12 
-0.95 


+0.20 
-1.17 


< -4.73 










UDS 


-4.91 


0, 


,31 


0.09 


+0.00 
-0.00 


+0.30 
-0.30 


+0.44 
-0.44 


< -5.21 










Total 


-4.29 


0, 


,13 


0.06 


+0.00 
-0.39 


+0.14 
-0.58 


+0.21 
-0.72 


< -5.33 







1.2 < z < 1.5 



2.0 - 2.1 



2.1 - 2.2 



2.2 - 2.3 



2.3 - 2.4 



2.4-2.5 



2.5 - 2.6 



2.6 - 2.7 



Cosmos 


— 4.11 


0.19 


0.08 


UDS 


-3.51 


0.06 


0.05 


Total 


-3.60 


0.05 


0.05 


Cosmos 


— 0.41 


u.uy 


U.UO 


UDS 


—3.12 


0.04 


0.04 


Total 


—3.18 


0.03 


0.04 


Cosmos 


-3.27 


0.07 


0.06 


UDS 


-3.02 


0.03 


0.04 


Total 


-3.07 


0.03 


0.04 


Cosmos 


-3.20 


0.07 


0.06 


UDS 


-3.03 


0.03 


0.04 


Total 


-3.07 


0.03 


0.04 


Cosmos 


-3.41 


0.09 


0.06 


UDS 


-3.24 


0.04 


0.04 


Total 


-3.27 


0.04 


0.04 


Cosmos 


-3.76 


0.13 


0.07 


UDS 


-3.94 


0.09 


0.06 


Total 


-3.89 


0.08 


0.05 


Cosmos 


-4.20 


0.22 


0.09 


UDS 


-4.99 


0.31 


0.09 


Total 


-4.63 


0.18 


0.06 



+0, 
-0, 


,00 
.10 


+0, 
-0, 


,34 
,10 


+0, 
-0 


,40 
,25 


—2.95 





,05 


0, 


,05 


+0, 
-0, 


,02 
,01 


+0, 
-0, 


,02 
,03 


+0.08 
-0.08 


+0, 
-0 


,04 
,01 


+0, 
-0 


,22 
,19 


+0, 
-0 


,23 
,20 


— Z. 1 1 


U 


no 
,Uz 


U, 


,Uo 


+0, 
-0, 


,00 
,07 


+0, 

-0, 


,03 
,05 


+0.05 
-0.09 


+0, 
-0, 


,04 
,01 


+0, 
-0, 


,21 
,18 


+0, 
-0 


,23 
,20 


—2.76 


0, 


,02 


0, 


,03 


+0, 
-0, 


,00 
,06 


+0, 
-0, 


,03 
,05 


+0.05 
-0.09 


+0 
-0 


,00 
,02 


+0 
-0, 


,06 
,36 


+0 
-0, 


,12 
,37 


—3.04 





,06 


0, 


,05 


+0, 
-0, 


,08 
,00 


+0, 

-0, 


12 

,16 


+0.16 
-0.18 


+0. 


!io 


+0 
-0, 


It 


+0 
-0 


,21 


-2.81 


0, 


,03 


0, 


,04 


+0, 

-0, 


,00 
,07 


+0, 

-0, 


,05 
,05 


-OXil 


+0. 


,00 
.09 


+0 

-o' 


OS 
.16 


-0 


,09 
^20 


-2.86 


0, 


,02 


0, 


,03 


+0, 
-0, 


,00 
,06 


+0, 

-0, 


,05 
,05 


+0.07 
-0.09 


+0, 


,01 
,13 


+0, 
-0 


,05 
,07 


+0, 
-0, 


,11 
,18 


-3.29 


0, 


,08 


0, 


,06 


+0, 
-0, 


,00 
,07 


+0, 
-0, 


,08 
,16 


+0.13 
-0.20 


+0, 


,00 
,05 


+0, 
-0, 


,01 
,02 


+0, 
-0, 


,05 
,07 


-3.00 


0, 


,03 


0, 


,04 


+0, 
-0, 


,06 
,09 


+0, 

-0, 


14 
,24 


+0.16 
-0.26 


+0, 


,00 
,04 


+0, 
-0, 


,01 
,02 


+0, 
-0 


,05 
,07 


-3.06 


0, 


,03 


0, 


,04 


+0, 

-0, 


,06 
,08 


+0, 

-0, 


13 
,23 


+0.15 
-0.24 


+0, 
-0, 


,08 
,00 


+0, 
-0, 


,00 
,03 


+0, 
-0 


,12 
,10 


-3.63 


0, 


,11 


0, 


,07 


+0, 

-0, 


12 
,00 


+0, 

-0, 


19 

,33 


+0.26 
-0.36 


+0, 


,00 
,06 


+0, 
-0 


,03 
,07 


+0, 
-0 


,06 
,11 


-3.38 


0, 


,05 


0, 


,04 


+0, 
-0, 


,03 
,05 


+0, 
-0, 


14 

,18 


+0.15 
-0.19 


+0, 


,01 
,06 


+0, 
-0, 


,02 
,07 


+0, 
-0 


,05 
,10 


-3.43 


0, 


,04 


0, 


,04 


+0, 
-0, 


,03 
,05 


+0, 
-0, 


13 
,17 


+0.15 

-0.18 


+0, 
-0 


,23 
,06 


+0, 
-0, 


,17 
,12 


+0, 
-0 


,30 
,17 


-4.51 


0, 


,31 


0, 


,10 


+0, 
-0, 


,48 
,00 


+0, 
-0, 


,54 
,00 


+0.79 
-0.32 


+0, 


,00 
,09 


+0, 
-0, 


,13 
,26 


+0, 
-0, 


,15 
,29 


-3.99 


0, 


,10 


0, 


,06 


+0, 

-0, 


11 

,05 


+0, 

-0, 


,43 
,60 


+0.46 
-0.61 


+0, 


,02 
,09 


+0, 
-0 


,12 
,25 


+0, 
-0 


,14 
,27 


-4.07 


0, 


,09 


0, 


,06 


+0, 
-0, 


11 

,04 


+0, 
-0, 


42 

,58 


+0.45 
-0.59 


' +0, 
-0 


,16 
,00 


+0, 
-0, 


,24 
,34 


+0, 
-0 


,32 
,37 


-4.81 


0, 


,43 


0, 


,11 


+0, 

-0, 


,00 
,00 


+0.30 
-99.00 


+0.54 
-99.00 


+0, 


,00 
,02 


+0, 
-0, 


,44 
,44 


+0, 
-0 


,46 
,45 


-4.99 


0, 


,31 


0, 


,09 


+0.00 
-99.00 


+0.40 
-0.00 


+0.51 
-99.00 


+0, 


,02 
,02 


+0, 
-0, 


,38 
,38 


+0, 
-0 


,39 
,39 


-4.93 


0, 


,25 


0, 


,08 


+0.00 
-99.00 


+0.34 
-99.00 


+0.44 
-99.00 


+0, 


,10 
,60 


+0, 
-0, 


,10 
,30 


+0, 
-0, 


,27 
,71 


< -4.81 




















tl 


,00 
,30 


+0, 
-0, 


,60 
,30 


+0, 
-0, 


,68 
,53 


-4.99 


0, 


,31 


0, 


,09 


+0.00 
-0.30 


+0.00 
-99.00 


+0.32 
-99.00 


+0, 
-0, 


,04 
,30 


+0, 
-0, 


,36 
,22 


+0, 
-0 


,42 
,43 


-5.11 


0, 


,31 


0, 


,09 


+0.00 
-0.30 


+0.00 
-99.00 


+0.32 
-99.00 



D. VELOCITY DISPERSION FUNCTION FOR STAR-FORMING GALAXIES - INCLUDING GAS MASSES 

As discussed in §2.2, increased specific star formation rates at z > 1 imply that galaxies at high redshifts must have 
elevated gas surface densities and therefore it may be important to include gas masses in addition to in calculating 
inferred velocity dispersions. For each galaxy, we calculate the inferred velocity dispersion using Mbaryon, as defined 
in §2.2 and recalculate the VDF. This has very little effect on the VDF for quiescent galaxies, but it does increase the 
velocity dispersions of star- forming galaxies, slightly shifting the number density distribution (by at most ~ 10%). In 
Figure [TT] we show the VDF for star-forming galaxies with (solid line) and without the inclusion of gas. These shifts 
are within the quoted error bars for the VDF and given uncertainties in calculated sSFRs, we ignore this effect for the 
analysis in this Paper. We expect the effect of gas to be increasingly important at higher redshifts than probed by 
this study. 
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Fig. 11. — Velocity dispersion function for star-forming galaxies calculated without (dashed lines and shaded polygons) and with (solid 
lines) inclusion of gas mass. 
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